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NONINTERSECTING BROWNIAN EXCURSIONS 

By Craig A. Tracy 1 and Harold Widom 2 

University of California-Davis and University of California- Santa Cruz 

We consider the process of n Brownian excursions conditioned 
to be nonintersecting. We show the distribution functions for the 
top curve and the bottom curve are equal to Fredholm determinants 
whose kernel we give explicitly. In the simplest case, these determi- 
nants are expressible in terms of Painleve V functions. We prove that 
as n — > oo, the distributional limit of the bottom curve is the Bessel 
process with parameter 1/2. (This is the Bessel process associated 
with Dyson's Brownian motion.) We apply these results to study the 
expected area under the bottom and top curves. 

1. Introduction. A Brownian excursion (BE) (AT(r)) Tg [o i i] is a Brownian 
path conditioned to be zero at times r = and r = 1 and to remain positive 
for < t < 1. This heuristic definition of this Markov process requires elabo- 
ration [11, 21] since such paths have zero Wiener measure. Alternatively, BE 
can be defined by scaling one-dimensional simple random walk conditioned 
to stay positive and conditioned to start and to end at the origin — a process 
known as Bernoulli excursion [22]. 

In this paper, we consider the process of n Brownian excursions condi- 
tioned to be nonintersecting, a model first introduced by Gillet [14] and Ka- 
tori et al. [16]. (Actually, the latter authors consider the closely related model 
of meanders instead of excursions.) Gillet proved that this n-nonintersecting 
BE is the distributional limit of suitably rescaled n Bernoulli excursions 
conditioned not to intersect. This Bernoulli exclusion model is known in 
the physics literature as vicious walkers with a wall and was first proposed 
by Fisher [9] as a statistical physics model of wetting and melting. Using 
the theory of symmetric functions, Krattenthaler, Guttmann and Viennot 
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[17] give explicit enumeration formulae for possible configurations of these 
walkers. 

If the BE paths are labeled Xi(t) < ■ ■ ■ < X n (r), < r < 1, then our main 
focus is on the lowest path, X±(t) and the highest path, X n (r). In particular, 
we show for each fixed positive integer m, with < r± < ■ • ■ < r m < 1, that 

P(x 1 (n)>x 1 ,...,x 1 (T m )>x m ) 

and 

P(X n (T 1 )<xi,...,X n (r m )<x m ), Xi £R + , 

equal Fredholm determinants of the form det(I — K be xj)i where K BE , the 
extended BE kernel, is an m x m matrix kernel and \J denotes multiplication 
by a diagonal matrix diag(xj fc ) whose entries are indicator functions. For 
the lowest path, J k = (0,x k ) and for the highest path, J k = (x k ,oo). 

The matrix kernel K BE (x, y) = (K®p(x, y))™ £=1 will have entries given by 
differences, 

K kF( x 'V) =H k e(x,y) - E M (x,y). 

Here, 

(1.1) E ke (x,y)-^ otherwise, 

where P-(x,y,r) is the transition probability for Brownian excursion, 

P _ ( Xj y> T ) = [L e -^W)/(2r) ginh ( *lA 
V 7TT \ r / 

The summand H k g(x,y) will have different representations, one of which is 
H k g(x,y) 



xV(2(l-r fe ))-3/V(2r,) 



T £ (l - Tfc) 

T"fe(l-^) 



xE 

i=o 



P2j+1 7= y= v P2j+1 



V^I-t^^ V v/2r,(l - r,) 

where the are the normalized Hermite polynomials. Although rather com- 
plicated, it has the same ingredients as the kernel for the Gaussian unitary 
ensemble (GUE). 

Extended kernels such as these were first introduced by Eynard and Mehta 
[8] in the context of a certain Hermitian matrix model. In Section 2, we give 
a self-contained account of a class of extended kernels. Though the results 
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appear in one form or another in the literature [8, 13, 26], we hope the 
reader will find this presentation beneficial. In Section 3, we use this theory 
of extended kernels to derive our expressions for the entries of K BE (x,y). 
Using the (scalar) m = 1 kernel, we show in Section 4 that for r £ (0, 1), 



P(X x (r) > s v /2r(l - r)) = exp (- J* ^dt 
where r = r(s) is the solution to 

(1.2) s 2 {r"f = 4(r') 2 (sr - r) + (sr' - rf - {An + l)(sr' - r)r' + \{r'f 
satisfying the asymptotic condition 

r(s) = r s 3/2 + O(s 5 / 2 ), s^O, 

where 

_1 1_ (2n\ An(2n + 1) 



A similar result is obtained for ¥{X n {r) < s^2t(1 — r)). 

The differential equations are the so-called a form of Painleve V equa- 
tions and the Fredholm determinants in these cases, where J= (0,x) or 
J = (x,oo), are Painleve V r-functions [12]. 

To obtain (1.2), we use results from [25] to obtain a system of differen- 
tial equations whose solutions determine r(s) and then reduce these to the 
second order equation for r(s). 

In Section 5, we consider the limiting process, as n — > oo, of a suitably 
rescaled bottom path X\(t) near a fixed time r £ (0, 1) and show that the 
distributional limit is ,61/2 ( r )) the Bessel process of [26] with parameter a = 
1/2. (This is different from the three-dimensional Bessel process appearing in 
the literature, which is Brownian motion conditioned to stay positive forever. 
When one rescales BE near r = 0, the three-dimensional Bessel process is 
the scaling limit.) This is a process whose finite-dimensional distribution 
functions are given in terms of an extended Bessel kernel K Bcs . Precisely, if 
-co < n < • • • < r m < 00, then 

(1.3) P(£i/ 2 (ti) > xi, . . .,B 1/2 (r m ) > x m ) = det(/ - K Bcs X j), 

where x.J — diag(xj fe ) with Jj. = (0, x^). For a = 1/2, the entries of the ma- 
trix kernel are 



/ e^x-^^smxtsmytdt, if k > 

00 

J e( Tk - T ^ t2 / 2 smxtsmytdt, if k < 
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Our result is that for fixed r € (0, 1), if, in the probability 

¥(X 1 (T 1 )>x 1 ,...,X 1 (r m )>x m ), 
we make the substitutions 



/r(l-r) r(l-r) 

then the resulting probability has the n — ► oo limit (1.3). This is proved by 
establishing trace norm convergence of the extended kernels. As far as we 
know, this is the first appearance of the Bessel process in the literature as 
the distributional limit of a finite-n Brownian motion process. 

We expect the top path, X n (r), again suitably rescaled, to have the Airy 
process [13, 20] as its distributional limit. The reason is that the top path, 
as n — > oo, should not feel the presence of the "wall" at x = and the Airy 
process is known to be the distributional limit of nonintersecting Brownian 
bridges. The bottom path does feel the presence of the wall and in the end 
this is why the Bessel process appears. We shall not derive the scaling limit 
at the top because it is the bottom scaling that gives something new, and is 
thus more interesting. 

For BE, a related random variable is the area A under the BE curve. 
(See [10, 18] and references therein.) Both the distribution of A and the 
moments of A are known (the moments are known recursively) and, perhaps 
surprisingly, numerous applications have been found [10, 19] outside their 
original setting. It is natural, therefore, to introduce A n ^L (A n ^n) equal to 
the area under the lowest (highest) curve in n-nonintersecting BE. From our 
limit theorem specialized to the m = 1 case, we deduce in Section 5.4 the 
asymptotics of the first moments E(^4 ni £) and E(-A n> #). In particular, 

E^Vl) ~ -^=, n^oo, 

where 

C L = -^J Q det(J - K Bcs X{ o,s)) ds ~ 0.682808. 

Here, K Bes is the m = 1 Bessel kernel, 

r d„, 2 f 1 , 1 /sin(x — y) sin(x + w)\ 

K (x, y) = — I smxtsmytdt= - ^ — i — . 

tt J 7r V x -y x + y J 
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2. Extended kernels. 

2.1. Initial and final points all distinct. Although only a special case of 
the following derivation will be used in our discussion of nonintersecting BE, 
it might be useful for the future, and is no more difficult, to consider a rather 
general setting. 

Suppose that {Xj(r)}™ =1 is a family of stationary Markov processes with 
continuous paths, with common transition probability density P(x,y,r), 
conditioned to be nonintersecting, to begin at a\, . . . ,a n at time r = and 
to end at b\ , . . . , b n at time r = 1 . An extended kernel is a matrix kernel 
K(x,y) = (K M (x,y)) k i =lj ,,, jin depending on < t\ < ■ ■ ■ < r m < 1 with the 
following property. Given functions f k , the expected value of 

m n 

nn^+AW 7 *))) 

k=\i=l 

is equal to det(I + K /), where / denotes multiplication by diag(/&). In the 
special case where f k = ~XJ k j this is the probability that for k = 1, . . . , m, no 
path passes through the set at time Tfe. For BE, this gives the statement 
made in the third paragraph of the Introduction, once we have computed 
the extended kernel for it. 

Here is how an extended kernel is obtained. Define the nxn matrix A by 

Aij = P(a,i,bj,l), 

define the m x m matrix kernel E(x,y) by 

(2.1) E ke (x,y)-^ otherwise, 

and if A is invertible, define the m x m matrix kernel H (x, y) by 

n 

(2.2) Hkt(x,y) = Pfab^l-TkXA-^jiPfaoiiTt). 

We shall show that H — E is an extended kernel in the above sense. 

From Karlin and McGregor [15] and the Markov property, we know the fol- 
lowing. For paths starting at xoj = a\ at time r = 0, the probability density 
that at times (k = 1, . . . , m + 1) the paths are at points Xki (i = 1, • • ■ , n) 
is a constant times 

m 

Y[ det(P(x kji ,x k+1>j ,T k+1 - Tk))i<i,j< n - 

k=0 

In the above, we set x m+ ij = bj,r m+ i = 1. The expected value in question 
is obtained by multiplying 

m n 

nn( i+ /*(^)) 

fc=li=l 
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by the probability density and integrating over all of the x k i- We apply the 
general identity 

det(ipj(x k ))™ k=1 ■ det(^j(x fc ))j ifc =i^( x i) ' ' ' dfj,(x n ) 
n!det( j tpj(x)ip k (x) d/j,(x) 

■ j,k=i 

successively to the variables xu, . . . ,x m i and find that the integral equals a 
constant times the determinant of the matrix with i,j entry 

„ „ m— 1 

/••• / P{Oi,XX,r{) Y[ P(Xk,X k+ l,T k+ i-T k )P(x m ,bj,l-T m ) 
k=l 

m 

x II ( X + fk{x k )) dx± ■ ■ ■ dx m . 
k=l 

If we set / = 0, the matrix is A by the semigroup property of P(x,y,r) 
and the expected value is 1, so the normalizing constant is (det^) -1 . Hence, 
if we replace P(ai,x,r) by 

Q(ai,x,r) = J2(A~ 1 ) ij P(a j ,x,T), 



then the expected value equals the determinant of the matrix with i,j entry 

„ „ m— 1 

/ •"• / Q( a U x l>Tl) II P ( X k,X k +l,T k+1 -T k )P(x m ,bj,l-T m ) 
k=l 

m 

x II ( X + fk{xk)) dxf- dx m , 

k=l 

which equals the determinant of / plus the matrix with i,j entry 

„ „ m— l 

/ ■•• Q(a i ,x l ,Ti) Y[ P{x k ,x k+ i,r k+1 -T k )P(x m ,bj,l ~T m ) 



k=l 



\[(l + fk(x k ))-l 



k=l 



dx\ ■ ■ ■ dx m . 



The bracketed expression may be written as a sum of products, 

r>l fci<---<fc r 

Correspondingly, the integral is a sum of integrals. By integrating with re- 
spect to all x k with k ^ kx, . . . ,k r , we see that the corresponding integral is 
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equal to 

/ • • • / Q(a,i,x kl ,T kl )f kl (x kl ) 




x P(x kl ,x k2 ,T k2 - r kl )f k2 (x k2 ) ■ ■ ■ P(x kr ,bj,l - T^) dx kl ■■■ dx kr . 

If we let U kj £ be the operator with kernel U k ^(x, y) = P(x, y, ti — T k )fi(y), 
then the above may be written as the single integral 



(If r = 1, we interpret the operator product to be the identity.) Replacing 
the index k\ by k and changing notation, we see that the sum of all of these 
equals 



^2Q(a,i,x,T k )f k (x)l U k}kl U kuk2 - ■■U kr _ ukr P(x,b j ,l-T kr ) dx, 



where the inner sum runs over all k r > ■ ■ ■ > k% > k. (Now, if r = 0, we 
interpret the operator product to be the identity.) 

Write Pj for the vector function {P(x,bj, 1 — Tfc)}^!^ and U for the op- 
erator acting on vector functions with matrix kernel (U k ^(x, y))^/^ , where 

we set U k /(x,y) = if k > I. Using the fact that (/ — U)^ 1 = J2 r >o U r , we 
see that the expression in large parentheses in the integral equals the kth. 
component of the vector function (J — U)~ 1 Pj. Denote it by T(k,x;j), and 
similarly denote the first factor in the integral by S(i;k,x). We may think 
of T(k,x;j) as the kernel of an operator T taking sequences into vector 
functions, 



while S(i;k,x) is the kernel of an operator taking vector functions to se- 
quences, 



The integral is the kernel of the product ST, evaluated at 

We use the general fact that det(/ + ST) = det(I + TS). The operator TS, 
which takes vector functions to vector functions, has matrix kernel with k,l 
entry 




P(x,bj,l - r kr ))dx. 







n 



J2T(k,x;j)S(j;£,y). 



.7 = 1 
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Since U is strictly upper-triangular, we have det(J — U) = 1, so 

det(J + TS) = det(J -U + (I- U)TS). 

If we recall that T(k,x;j) is the kth component of (I — U)~ l applied to the 
vector function {P(x, bj, 1 — r k )} and recall the definition of S, we see that 
(/ — U)TS has k,l entry 

n 

J2P(x,bj,l-T k )Q{aj,y,T £ )f e (y) 

3=1 

n 

= XI ^A^ 1 - r fc)(^~% p (y> a *'' ^ ^).My)• 
£J=l 

This equals Hki(x,y)fi(y). Since £7 has k,£ entry E k i(x,y)fi(y), we have 
shown that the expected value in question equals det(I + (H — E)f). Thus 
H — E is an extended kernel. 



2.2. Initial and final points at zero. In the preceding derivation, we as- 
sumed that A was invertible. If any two are equal or any two bj are equal, 
this will certainly not hold. There are analogous, although more complicated, 
formulas if several of the Oj or bj are equal. We consider here the simplest 
case where they are all equal to zero. We take (2.2) for the case when they 
are all different (when A is expected to be invertible) and compute the limit 
as all of the and bj tend to zero. 

There is a matrix function Dq = Dq(c\, . . . ,c n ) such that for any smooth 
function /, 



/ /(0) \ 



lim Dt) 



f(02 



\f(Cn)J 



\ f(n-l) (0)7 



Here, lim Ci ^ c is short for a certain sequence of limiting operations. The 
matrix D' (the transpose of Dq) acts in a similar way on the right on row 
vectors. 

Writing D a for D$(ai, . . . , a n ) and for Do(bi, . . . , b n ), we use 

A- l = D' b (D a AD' b y l D a 

and then take limits as all ai,bj — ► 0. The matrix D a AD' b becomes the matrix 
B with i,j entry 

B ij = di- 1 di- 1 P(a i ,b j ,l)\ ai=bj=0 , 
the row vector {P(x, bj, 1 — r^)} becomes the vector with jth component 

1 - T k ) = dl^Pix, bj,l- T k )\ bj=0 
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and the column vector {P(ai,y,Te)} becomes the vector with ith component 

P?(y,n) = di-: 1 P(a j ,y,n)\ ai=0 . 

Then if the matrix B is invertible, we obtain for the iif-summand of the 
extended kernel, in this case, 

n 

(2.3) H M (x,y)= Pjfal-TkHB-^jilffaTi). 

*.i=i 

One would like B to be as simple as possible. Something one could have 
done before passing to the limit would be to take two sequences {at} and 
{f3j} of nonzero numbers and in (2.2), to replace A by the matrix (otiAijPj), 
replace the function P(x,bj,l — 7*) by P(x,bj,l — Tk)f3j and replace the 
function P(ai,y,Ti) by aiP(ai,y,re). The kernel does not change. 

2.2.1. An example. We apply (2.3) to the case of n nonintersecting Brow- 
nian bridges. The transition probability for Brownian motion is 

P(x,y,T) = ^e-^ 2 ^ 

and we choose at = e a ^ 2 and (3j = e b ^ 2 . Then the modified A equals (e aib i j 
\/2~7r), the modified B equals diag((z — l)!)/\/2~7r, the modified B -1 equals 
V2tt X diag(l/(t- 1)!) and 



-x 2 /(2(l-r fc )) 



^(l-Tfc) 



(j'-l)/2 



l 



2(1-7*; 

_ p -?/ 2 /(2ri) 



Hj-i 

1-T£ 



X 



2r £ 



V2r fc (l-r fc ) 

(i-l)/2 



a; 



\/27rT£ 

where H k (x) are the Hermite polynomials satisfying 



i-1 



V2r £ (l-r £ ) 



fc>0 



In terms of the harmonic oscillator wave functions, 



(2.4) 
we have 



H u (x,y)= ... 1 . e (V 2 -^ 



^2(1-7*)^ 



n-l 
j=0 



Tfc(l -Tf) 
7jf(l - T k ) 
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where, for notational simplicity, we write Xk = xj \J2t]~(\ — t\) and = 
yj \I2ti(\ — Tg). The matrix kernel H(x,y) — E(x,y) is called the extended 
Hermite kernel. For m = 1, Hn(x,y) (with r± = r) reduces to the Hermite 
kernel K^ E (X, Y). [We can neglect the two exponential factors since they 
do not change the value of the determinant. The factor \j \J2t{\ — r) also 
disappears when we multiply by the differential dY .] [In the usual defini- 
tion of the extended Hermite kernel, the times 77% € (0, 1) are replaced by 
times ffc G (0,oo) defined by Tfc/(1 — r^) = e 2l ~ k so that H — E has k,£ entry 



K^ E (X,Y;f k -f e ), where 



(2.5) 



K^ E (x,y;f) 



n-l 
j=0 



if f > 0, 



^2e JT ip j (x)<p j (y), iff<0. 



The representation for f < follows from the Mehler formula ([7], Section 
10.13 (22)).] 

2.3. General case. In the recent literature, one finds investigations where 
not all starting and ending points are the same. In [2, 27], there is one 
starting point and two ending points and in [6], there are two starting points 
and two ending points. 

The previous discussion can be extended to the case where there are 
R different and S different bj. The different a^, we label a±, . . . ,or and 
assume that a r occurs m r times. Similarly, the different bj, we label 61, . . . , 65 
and assume that b s occurs n s times. Of course, Y^ m r = J2 n s = n. It is 
convenient to use double indices ir and js where i = 0, . . . , m r — 1 and j = 
0, . . . ,n s — 1. We may think of A as consisting of m r x n s blocks A rs . 



Once again, there is a matrix function D c = D c (c\, 

( m \ 

f'(c) 



, c„ ) such that 



lim D r 



/(<*) 



\f(Cn)/ 

and similarly for the action of D' c on row vectors. We then write 

where D a and D b are now block diagonal matrices, 

D a = diag (D ar ) , D b = diag (D bs ) , 

substitute into (2.2) and take the appropriate limits. The matrix entry A{ r j t 
is replaced by 



Bir,j s — d a . r di P(ai r ,bj s , 1) 



-a r ,bj 



NONINTERSECTING BROWNIAN EXCURSIONS 



11 



while the row vector {P(x,bj s , 1 — r^)} becomes the vector with js compo- 
nent 

P jt fa. 1 - T k) = dt^Pfa b ja , 1 - T k )\ bjs=bs 
and the column vector {P(oj r , y, r^)} becomes the vector with ir component 

P-; (y, n) = di- l p{a ir , y , n)\ air=ar . 

Thus (2.3) is now replaced by 

H ke (x,y) = J2 P-KxA-r^B-^irP^iy,^). 

ir,js 

3. Nonintersecting Brownian excursions. 

3.1. Preliminary result. Let B T denote standard Brownian motion (BM), 
so that 

F X (B T G dy) = P(x, y, r) dy = _J =e -(*-») a /(2r) d y, r > 0, x, y G K, 

V 27TT 

and the subscript x on P indicates that the BM was started at x. Let 

H a = infjs : s > 0, B s = a} = hitting time of a. 
The basic result we use is [11, 21] 

(3.1) P X (B T G dy, H >t) = (P(x, y, r) - P(x, -y, r)) dy =: P. (x, y, r) dy. 

Equation (3.1) gives the probability that a Brownian path starting at x is in 
the neighborhood dy of y at time r while staying positive for all times < 
s < t. This process has continuous sample paths [11]; therefore, by Karlin- 
McGregor [15], we can construct an n nonintersecting path version and apply 
the formalism of Section 2. It is convenient to write 

(3.2) P_ (x, y, t) = \ /A e -(*W)/(2T) sinh 

V 7TT 

[The theory of Section 2 applies equally well to the n nonintersecting path 
version of reflected (elastic) BM, Y T := \ B T \. In this case, 

Vx(\B T \ G dy) = {P(x, y, r) + P(x, -y, r)) dy 

= ,/A e -(^ + , 2 )/(2.) cosh ^V 

V 7TT V T J 

We leave it as an exercise to construct the extended kernel.] 
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3.2. Extended kernel for Brownian excursion. Referring to (2.2) and the 
technique described after (2.3), we find that for distinct <n and bi, 

— n 

I J2 P-{xMA--n t )e^ / \A- 1 ) ij P-(y,a j ,T t )e^' 2 , 

where Aij = sinh(aj6j). The even order derivatives of sinhx all vanish at 
x = 0, so to take a^bj — > 0, we have to modify the procedure described in 
Section 2.2. (Otherwise, the matrix B so obtained would be noninvertible.) 
We now use the fact that there is a matrix function D c - such that for any 
smooth odd function /, 

/7(d) \ / f(0) \ 

f(02) /'"(0) 



lim D c 



\/(c„)/ 



V/(2n-i) (0) y 



Thus the B matrix is 

Ba = dg- 1 ^- 1 sinh(a i 6 i )k=6 J -=o = (2j - 1)% 
It follows that the limiting kernel has k,£ entry 



(3.3) 



, -re— 1 -i 
'7T x - 1 



H ki ( x ,y) = (2j + i)\ Pj ( x,1 ~ T k) p 3{y^n) 



where 

P j (x,T) = [dy +1 P^x 1 y,T)ey 2 /% = o 

= 1 c -^ 2 /(2r) 

X [^J+^-CC^^/WJ^+^/r _^'+l e -((l-T)/(2T)) 2 / 2 -x 2 //T 



ly=o 



2 /1-r 



7tt V 2r 



j+l/2 



V2r(l-r) 



e -^/(2r )jff2 . +l| 

Thus the ii-summand of the extended kernel can be written as 

-x 2 /(2(l-r k ))- y y(2r e ) 



H k e(x,y) 



n(i - Tk) 



(3.4) 



3=0 



Tk( 



n(i - r k ) 



P2j+1 



^2r k {l-T k ) 



P2j+1 



y/2Te(l-T t ) 
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Recall that p k are the normalized Hermite polynomials. It follows from the 
Christoffel-Darboux formula (see, e.g., [7], Section 10.3 (10)) for Hermite 
polynomials that 

n-l 

J2P2j+l(x)p2j+l(y) 
3=0 

U f p 2n (x)p2n-l(y) ~ P2n-\ {x)p2n (y) 



2 I x-y 

P2n(x)p2n-l(y) + P2n-l(x)p2n(y) 



+ 



x + y 



Setting Tfc = r, X = x/^2t(1 — r) and Y = yj \J2t{\ — r), we see that the 
diagonal elements are 



Hkk{x > y) = ^2W^ e 



-x 2 /(2(l-r))-j / 2 /(2r) 



(3.5) 



P2n{X)p2n-l{Y) - P2n-\ (X)p 2n QQ 

X-Y 

P2n(X) P 2n-l(Y) + p 2 n-l (X)l>2n (Y) 



+ 



e (l/2-r)X 2 e -(l/2-r)y 2 



X + Y 

{K^ E (X,Y)-K^ E (X,-Y)}. 



V2r(l-r) 

More generally (recall the notation of Section 2.2.1), 

(l/2-n,)X*-W2-T t )Y? 

H M (x, y) - E ke (x, y) = {^2 G n UE (*fc, Y t ; f k - f e ) 

(3.6) 

- K 2n E (^k, -Ye;f k - f e )}, 

where K^ E (x, y; f ) is defined as in (2.5). (To establish this identity for 
k < £, one uses the (slightly modified) Mehler formula, 

00 

J2 9 2j+ V2J+1 (x)f2j+i (y) 
3=0 

e -(i+ 9 2 )/(2(i- g 2 ))(x 2 +j / 2 ) sinh ^ 2 ^y 



" V^-q 2 ) U-<zV 

with x,y -^X k ,Y e and q = [r fc (l - n)Jn(l - n] 1 / 2 .) 
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3.2.1. Density of paths. The density of paths at time r is given by 
p n (x,r) = Hn(x,x) 

(3-7) = * {K^ E (X, X) - K^ E (X, -X)} 

V2r(l-r) 

cy n—1 
V2t(1 - r) J=0 
and satisfies the normalization 



p n (x, r) da; = n. 

For n = 1, p n reduces to the well-known density for Brownian excursion. 

4. Distribution of lowest and highest paths at time r. 

4.1. Fredholm determinant representations. Recall that we assume the 
nonintersecting BE paths to be labeled so that for all < r < 1, 

X n (r)>--->Xi(r)>0. 

We are interested in the distribution of the lowest path, X\ and the highest 
path, X n . From Section 2, we see that 

P(Xi(r) > x)=F(X 1 (t) >x,.. .,X n (r) >x)= det(7 - K XJl ), 

(4.1) 

Ji = (0,s), 

and 

P(X n (r) < x) = P(Xi(r) < x, . . . ,X n (r) < x) = det(7 - 7T X J 2 ), 

(4.2) 

7 2 = (s,oo), 

where, in both cases, 

(4.3) K(x,y) = K^ E (x,y) - K^ E (x, -y) 

and s = xj \/2t(1 — r). [Recall that the relation between K and K BE is 

K(x,y) = ^2t(1 - r)7C BE ( v /2r(l - t)x^2t{\ - r)y), 

aside from the exponential factor appearing outside the curly brackets in 
(3.6). For the determinants, we may ignore this factor.] 
The operator K is given by 

n-l 
3=0 
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Table 1 

The expected areas of lowest and highest curves for 
n=l,...,9 



n E(A n , L ) E(A n , H ) 



1 


=0.626657... 
f (V2-1)VtF= 0.458859... 


y | = 0.626657 
§V?r = 1-107783. 


2 


3 


0.380211. .. 


1.479479. . . 


4 


0.331955. . . 


1.791454. . . 


5 


0.298441. . . 


2.065097. . . 


6 


0.273407. . . 


2.311582. . . 


7 


0.253784. . . 


2.537567. . . 


8 


0.237865. . . 


2.747386. . . 


9 


0.224612. . . 


2.944040. . . 



where ^j(x) = V2(f2j+i(x) and where u ® v for functions u and v denotes 
the operator / — ► u(v, f). It follows that 

(4.4) det(J - K XJ ) = det(S jk - (%,* fe ))^l > 

where the inner product (•, •) is taken over the set J. Thus (4.4) gives finite 
determinant representations of P(Xi(r) > x) and P(X n (r) < x). 

As an application of (4.4), consider the areas under the lowest curve and 
the highest curve, 

K,L-= I X x {r)dT, A n , H := [ X n {r)dT. 
Jo Jo 

For n = 1, Ai l — A\^h an d the distribution of this random variable is called 
the Airy distribution (not to be confused with the Airy process!), the mo- 
ments of which are known recursively (see [10] and references therein). The 
expected area under the lowest curve is 

E(A n , L )= E(X 1 (r))dT= / P(Xi(r) > x) dxdr 
Jo Jo Jo 

(4.5) = j ^2t{1 - t) dr J™ det(/ - K X{0 , s) ) ds 

(4-6) = ^J det (*i.*-(*i'**))i,fci) tfe ' 

where the inner product is over the interval (0, s). A similar formula holds for 
E(A nt i{). Table 1 gives some data on E(A„ : i) and E(A nj #). Higher moments 
are expressible in terms of integrals over the Fredholm determinants of the 
extended kernels. We have not investigated these higher moments. 
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4.2. Integrable differential equations. In this section, we show that the 
Fredholm determinant of the integral operator with kernel 

(4.7) K(x,y) XJ (y) = [K^ E (x,y) - K^ E (x,-y)]xj(y) 

can be expressed in terms of a solution to a certain integrable differential 
equation when J = (0, s) and J = (s, oo). To do this, we transform the kernel 
K{x,y) to the canonical form 

/ 48 n V{x)ip{y) - i/>(x)<p(y) 

x-y 

so that the theory of Fredholm determinants and integrable differential equa- 
tions, as developed in [25], can be directly applied. If, instead of K(x,y), we 
use the kernel 

Ko{x ' y)= 2(^)1/4 ' 

then 

(4.9) det(J - K X{0>S 2)) = det(J - K oXi o, s ))- 

By use of the Christoffel-Darboux formula, we see that Kq is of the form 

(4.8) , with 

(4.10) rW/Wv^), V(^)=n 1/ V 1 /V2n-i(v / ^), 

where the (fk are the harmonic oscillator wave functions (2.4). 

Using the well-known formulas for the action of multiplication by x and 
differentiation by x on the (pi~, we find 

(4.11) xip'{x) = (4 — \x)tp{x) + y/nxi/)(x), 

(4.12) xip'(x) = -^ip{x)-(\-\x)ip{x). 

These formulas will allow us to use the results of [25]. 

In the derivation of the next section, we write J for (0, s) and denote by 
R(x,y) the resolvent kernel of Kqxj, the kernel of (/ — Kqxj)~ 1 Kq\J- Of 
course, this also depends on the parameter s. A basic fact is that 

(4.13) ^logdet(I -K oXJ ) = -R(s,s). 

Once R(s, s) is known, integration and exponentiation give the determinant. 
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4.2.1. Reduction to Painleve V, J = (0,s). Following [25], we define the 
functions of x (which also depend on s since J does) 

Q = (I-K XJ)- 1 >P, P=(I-K oX j)- 1 *P 

and then the quantities (which depend on s only) 

q = Q(s), p = P(s), r = sR(s,s), 

u = (Q,<p), v = (Q,ifj), w(s) = (P,ip), 

where the inner products are taken over J. 

The definitions are not important here. What is important are the rela- 
tions among the various quantities, as established in [25]. Of course, we are 
interested in r. This function and its derivative are given in terms of the 
other functions by 

r(s) = Cw + \/n)q 2 + (i — s — 2y/nw)qp 

(4.14) 

+ (—14 + \/ns + 2y/nv)p 2 , 

(4.15) ^ = -qp+ yfnp 2 . 
ds 

((4.14) is the specialization of (2.27) of [25] with additional insertions de- 
scribed in (2.32); (4.15) is the specialization of (2.28) of [25] with additional 
insertions described in (2.35). In both cases, we use the formulas (4.11) and 
(4.11).) 

Specializing the differential equations of [25] to the case at hand, we find 
that 

(4.16) s — = (\ - \s - \/nw)q + (—it + 2^fnv + \/ns)p, 

(4.17) s^ = —(w + \/n)q—(\ — \s — \/nw)p, 
. , , du 9 dv dw 9 

(4.18) -t = t, t=pq> ~r=r- 

ds ds ds 

((4.16) and (4.17) are (2.25) and (2.26) of [25], resp., with additional inser- 
tions described in (2.31). We also used recursion relations (2.12) and (2.13) 
of [25], evaluated at x = s. Equations (4.18) are "universal equations" given 
by (2.16)-(2.18) of [25].) 

This is a system of five equations in five unknowns. One has the first 
integrals, discovered by a computer search, 

(4.19) sp(q — y/np) + \/n(u + \w) — (1 + 2n)v + w{u — 2y/nv + nw) = 0, 

(\/ns — u + 2\fnv)p 2 + (y/n + w)q 2 

(4.20) 

+ (3 — s — 2\fnvS)pq + v — yfnw = 0. 
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(A first integral for a system of differential equations is a relation stating 
that a particular expression in the variables is constant when the system is 
satisfied. Often, as here, the constant can be determined by using the initial 
conditions.) 

To verify either of these, we differentiate the left-hand side using (4.16)- 
(4.18) and find that the derivative equals zero. Thus the left-hand side is a 
constant. That the constant is zero follows from the easily verified fact that 
all quantities are o(l) as s — > 0. 

In addition to (4.14) and (4.15), we compute r" using (4.16)-(4.18). These, 
together with the first integrals (4.19) and (4.20), give us five equations. 
Remarkably, from these, we can eliminate q, p, u, v and w and we obtain 
the second order equation (1.2) for r. (For the elimination, we used a Grobner 
basis algorithm as implemented in Mathematica. A Grobner basis of an 
ideal in a polynomial ring is a particular basis that is well suited for the 
elimination of variables. The reader is referred to the book by Cox, Little and 
O'Shea [5] for both the theory of Grobner bases and many examples of their 
application to computational commutative algebra. Remarkably, although 
one would expect to need six equations to eliminate the five unknowns, we 
are in a nongeneric situation and these five equations suffice. This would not 
have been discovered without the use of the Grobner basis. In our original 
derivation, we used only the first integral (4.19) and obtained a third order 
equation in r. Cosgrove [4] knew how to reduce this to the second order 
equation (1.2). This then led us to discover (4.20).) The solution we want 
satisfies the asymptotic condition 

(4.21) r(s) = r 

s 3/2 + q (s 5/ 2) 

as s — ► 0, 

where 

(4.22) 1 1 /2nA4n(2„ + l) 



This condition follows from the fact that R(s, s) ~ Kq(s, s) as s — > 0. (Power 
counting suggests that cs 1//2 is the leading term, but a computation shows 
that c is zero. The s 3//2 term is the leading one. Higher terms can be deter- 
mined from the differential equation.) 

4.2.2. Painleve V representation o/P(Xi(r) > x). From (4.13) and (4.9), 
we conclude that 

(4.23) P(Xi(r) > s^J2t{1 - t)) = exp(- jT ^dt 

where r is the solution to (1.2) satisfying boundary condition (4.21). Using 
(1.2), it is routine to compute the small s expansion of r(s) and hence that 
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of (4.23), 

P(Jf!(r) > s^M^7)) = 1 - \r s* + |(4n + l)r oS 5 - ^+^ + 9 ^ 

(4n + l)(32n 2 + 16n + 15) „ 

+ — -rns 

8505 ° 

128(2n + 3)(n-l) 2 10 » 

+ 275625 r ° S +U{S h 

where vq is defined as in (4.22). 

We remark that the number of BE curves, n, appears only as a param- 
eter in the differential equation (1.2) and the boundary condition (4.21). 
Although we have explicit solutions in terms of finite n determinants, ex- 
pansions such as the one above are more easily derived from the differential 
equation representation. 

4.2.3. Reduction to Painleve V, J = (s,oo). This differs from J = (0,s) 
in only minor ways. In this case, 

— log det(I - Koxj) = R(s, s), 

so that with r defined as before, 

nX n {r) < s^2t{1-t)) = exp(- jT ^ dt) . 

Specializing the differential equations of [25], we find for J = (s,oo), 



s = ( \ — \s — y/nw)q + (—u + 2yfnv + y/ns)p, 



du n dv dw 2 

ds ds ds 

The connection between r and the above variables remains the same as that 
given in (4.14). One has the two first integrals 

—sp(q — \/np) + y/n(u + \w) — (1 + 2n)v + w(u — 2\fnv + nw) = 0, 

(^/ns — u + 2^/nv)p 2 + (yfn + w)q 2 + (| — s — 2yfnw)pq — v + \fnvo = 0. 

Proceeding with the calculations as described above, we obtain 

s 2 (r") 2 = — 4(r') 2 (sr' — r) + (sr 1 — r) 2 — (4n + l)(sr' — r)r' + \{r') 2 . 
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5. Limit theorems as n — > oo. In this section, we obtain the scaling 
limit for the extended BE kernel at the bottom. We first explain what this 
means and then state the result. 

Given a kernel L(x, y), the kernel resulting from substitutions x — ► ip(x),y - 
<p(y) is the scaled kernel 

L(x,y) = (ip'(x)<p'(y)) 1/2 L(<p(x),(p(y)), 

which defines a unitarily equivalent operator. If there is an underlying pa- 
rameter and L has a limit Lq (pointwise or in some norm), then Lq is the 
scaling limit of L(x,y). If the convergence is in trace norm, then the limit 
of det(I — L) equals det(J — Lq). Trace norm convergence is assured if the 
underlying interval is bounded and if L(x,y) and its first partial derivatives 
converge uniformly. 

The extended Bessel kernel K B [24], Section 8, is the matrix kernel 
whose entries are given by 

1 e z2 ( Tk - Te)/2 <S> a (xz)$ a (yz) dz, if k > £, 



K%F S (x,y) 



/oo 



where $ a ( z ) = VzJa( z )- In particular, &i/2( z ) = \/2/7r sinz. 

We shall show the following. Take a fixed r G (0, 1) and make the substi- 
tutions 



and replacements 

(5.2) r fc ^r+ r(1 2 ~ T) r k . 

Then the extended BE kernel has the extended Bessel kernel with param- 
eter a = 1/2 as scaling limit as n — > oo, in trace norm over any bounded 
domain. In particular, we deduce the result stated in the introduction, that 
the finit-dimensional distribution functions for the BE process, when suit- 
ably normalized, converge to the corresponding distribution functions for 
the Bessel process. 



5.1. Integral representation for the extended kernel. The extended BE 
kernel K BE has the form (H^) — (Eke), where Ef-i is given by (1.1) and 
is given by (3.3). The limit theorem for the Hke will be derived from an 
integral representation of (3.3) which, as in [3, 27], uses two different integral 
representations for Pj. (We could use known integral representations for the 
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Hermite polynomials, but the method we use here is more general and also 
works in other cases.) 

The first is easy. With integration over a curve around zero, 



P_( x ,y,r) e y 2 / 2 = ^-Jp^x,t,r)e t2 / 2 - 



t-y 



so 



2vri 
2 



dt 

t 2jT2 



-x 2 / 



2r(2j + l)! / -ri 



-(l-)/(2r)t 2 sinh 



xt dt 



7TT 27TZ 

For the other, we take the Fourier transform on the y-variable. A compu- 
tation gives 



/oo 
e tay P-(x,y,T)e y /2 dy 
-oo 



= e - 2 /2(l-r) e -r/(2(l-r))^ sinh 



/.r.s 



^T^r~ ~ 1-r 

Taking inverse Fourier transforms and differentiating gives 
2 



Pi(x,r) 



3 x 2 /2(l-r) 



7T\/1 — T 

roc 



X 



/■°° e -r/(2(l-r))-» ginh J^ { _ is) 2j+1 ds 

J-oo 1-r 

/ZOO _ rr> c 

e r/(2(l-r))^ sinh ^ s 2,+l ds 
-ioo 1 r 



7r\/i — r 



,x a /2(l-r) 



These give 

H k t{x,y) 



V(l-^)(l-r/) 



-x 2 /(2(l-r fc )) +? y 2 /(2(l-r f )) 



7T 



i0 ° e - Tfc /(2(l-r fc ))t2 +T ,/(2(l-r,))s 2 sinh . 
-ioo 



x sinh ■ 



2/s 



2/i 



1 



s cis dt 

IF^t?' 



1-r, 

If we ensure that the s-contour passes to one side of the i-contour, then the 
summand 1 in the bracket contributes zero since the t-integral is zero. Thus 
we have 

1 



7rV(l-7*)(l-7*) 



-* 2 /(2(l-T fc )) +? , 2 /(2(l-T,)) 
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(5.3) X i / e 



-T k /(2(l-T k ))t 2 +T e /(2(l-T e )) S 2 



xt ys ( s\ 2n sdsdt 

x smh smn ■ 



1 — Tfc 1 — T£ \ t J s 2 — t 2 

5.2. Limit theorem for the bottom. If we ignore the denominator s 2 — t 2 
in (5.3), we have a product of two integrals for each of which we can apply 
the method of steepest descent. If we think of x and y as small and k = £ 
with Tfc = Te = r, then there are saddle points at ±i^/2n(l — tJJt for both 
integrals. So we assume that all are near a single r and make the changes 
of variable 



„ J** 1 -<■)., t ~f" ( \- T) t. 

This is what suggested the substitutions (5.1) and replacements (5.2). Ob- 
serve that in the scaling limit, t remains the same. 

For (5.3), we can ignore the outside factor involving x and y since remov- 
ing it does not change the determinant. An easy computation shows that 
after the various substitutions, the rest of Hki(x,y) becomes 

4 °° e n(s 2 -t 2 +21ogs-2logt)-(T k +0(n- 1 ))t 2 /2+(r e +0(n- 1 ))s 2 /2 



TT 2 + 0(n- 



xt ys sdsdt 

x smn — — t— smn - 



l + Oin^ 1 ) l + 0{n~ l ) s 2 -t 2 ' 

where the terms 0(n~ l ) are independent of the variables. 

The saddle points are now at ±z for both integrals. The steepest descent 
curve for the s-integral would be the imaginary axis and for the t-integral 
would be the lines Imt = ±1, the lower one described left to right, the upper 
one right to left. But in order to get to these steepest descent curves, we 
must pass the s-contour through part of the original t-contour. To be specific, 
assume that the t-contour was the unit circle described counterclockwise and 
that the s-contour passed to the right of the t-contour. Then if we replace the 
s-contour by the imaginary axis (after which we can replace the ^-contour by 
the lines Imt = ±1), we must add the contribution of the pole at s = t when 
t is in the right half-plane. Thus the change of contour results in having to 
add 

ZL^ f l e^+o(n-^t 2 /2 sinh ^_^ inh ^—dt. 

m + 0(n- l )J-i l + 0{n~ l ) l + 0{n~ l ) 

This is 



vr jo 



/ e {Tk ~ Te)t2/2 smxtsmytdt + 0(n~ 
Jo 
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uniformly for x and y bounded. 

It is easy to see that the remaining integral, the s-integral over the imag- 
inary axis and the t integral over Imi = ±1, is (9(n -1 / 2 ) uniformly for x,y 
bounded. In fact, outside neighborhoods of the points where s 2 = t 2 , we may 
estimate the double integral by a product of single integrals, each of which 
is 0(n -1 / 2 ), so the product is 0{n~ l ). There are four points where s 2 = t 2 , 
where both s and t are ±i or where one is i and the other —%. If we set up 
polar coordinates in neighborhoods based at any of these points, then the 
integrand, aside from the last factor, is 0(e~ Snr ) for some S > 0, while the 
last factor is 0(r _1 ) because s and t (or s and —t) are on mutually perpen- 
dicular lines. Therefore, because of the factor r coming from the element of 
area in polar coordinates, the integral is O(J °° e~ Snr dr) = 0{n~ 1 / 2 ). 

Thus Hki(x,y) has the scaling limit 

- / e iTk ~ Te)t /2 sinxts'mytdt 
vr Jo 

uniformly for x and y bounded. When k < £, we must subtract from this 

Eke(x,y) =P_(x,y,r£-r fc ), 

which equals the expression above, but with the integral taken over (0, oo). 
So the effect of the subtraction is to replace the expression by 

2 f°° 2 

/ e ( Tfc-T ^ /2 sinxtsmytdt. 

vr Jx 

This limiting matrix kernel is exactly the extended Bessel kernel K Bes (x,y) 
with a = 1/2. 

We have shown that the k,i entry of the scaled BE kernel converges to 
the corresponding entry of the extended Bessel kernel, uniformly for x,y 
bounded. A minor modification of the computation shows that the corre- 
sponding result holds for derivatives with respect to x and y of all orders and 
this is more than sufficient to give trace norm convergence on any bounded 
set. This concludes the proof of the result stated at the beginning of the 
section. 

5.3. Limit theorem for highest curve. As mentioned in the Introduction, 
we shall not derive the scaling limit at the top. But since below, we shall use 
the top scaling of the diagonal terms Hkk(x,y) to the Airy kernel, we point 
out that this is very easy. With the usual scaling X -> y/2n + X/(2 1 / 2 n 1 / 6 ), 
Y — > \^2n + Y/ (2 1 / 2 ?! 1 / 6 ), the summand in (3.5) with the denominator X + Y 
is easily seen to converge in trace norm to zero over any interval (s,oo), by 
the asymptotics of the Hermite polynomials. The rest of the kernel converges 
to the Airy kernel in trace norm, by known results. 
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5.4. Asymptotics of 1 &{A n ^) and E(A nj #). From our limit theorems fol- 
low the asymptotics of E,(A n ^) and E(^4 ni #). 

For the lowest curve, we use formula (4.5), which tells us that 

E{A n , L ) = — / det(J - Kxto,,)) ds, 



471 Jo WI-Kxm), 

where K is given by (4.3). It is given in terms of K BE with m = 1 and 
arbitrary t\ = r by 

K(x, y) = ^2t(1-t)K be (^2t(1-t)x, ^2t(1 - r)y) 

times exponential factors which we may ignore. Hence, the scaling we have 
derived shows that 

1 R ( y_ \ _^ k V*b/ x y ) = -( sin ^ ~ y "> sin ^ x + y "> 



2^fn \2- v /ra' 2^/n J ' ir \ x — y x + y 

in trace norm on bounded sets. As in Section 4.2, instead of K Bcs , we use 
defined by 



2xV4yl/4 



It is convenient to use K® cs since we can make direct use of the differential 
equations for it derived in [24]. 
By a change of variable, we have 



E(A„, L ) = -^= r°det(I-Kx 
8v2n Jo 



(0,s/2^/E)) ds - 



Therefore, if we can take the limit under the integral sign (which we shall 
show below), then this implies that 

(5.4) E(A, L )~^=, n^oo, 



where 

After expressing this in terms of K Bes , it becomes the formula quoted in the 
Introduction. 

Solving numerically the differential equation satisfied by the logarithmic 
derivative of det(7 — if Bes X(o,;c)) [24], followed by a numerical integration, 
gives 

r°° l 

/ -^det(I-K^ s X{0 , x) )dx^ 4.917948 
Jo \X 
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and hence cl — 0.682808. From Table 1, we compute y/nK(A n ^) for n = 
5,. ..,9: 

0.667334, 0.669708, 0.671449, 0.672784, 0.673838. 
To justify taking the limit under the integral sign in 

POO 

det(J - K X(o,s/2^)) ds > 
it is enough to show two things: 

/•OO 

(5.5) lim / det(I - Kx {0)S/V E))ds = 0, 

n 00 J v^nlogn 

det(I - Kx(o,s/2^)) 
(5-6) _ 

= 0(e s ) uniformly for some 5 when s < y/nlogn. 

Combining (5.5) with an application of the dominated convergence theorem 
using (5.6) justifies taking the limit under the integral sign. This completes 
the proof of (5.4). 

The eigenvalues of -f^X(o,s/2 v / H) ue m the interval [0, 1] because K is a 
projection operator. For (5.5), we use the fact that the largest eigenvalue of 

K X(0,s/2y/n) is at least that of 2i Pl ® PlX(0,s/y/E)i so 



rs/2^/n poo 

det(I-2fx(o,«/Vn)) < 1-2 / ifi{xfdx= fi{x) 2 dx = 0(e 

J0 Jsl2jn 



s 2 /4n\ 

s/2^n 



Integrating this over (-^/nlogn, oo) gives (5.5). 

For (5.6), we use the fact that since the spectrum of Kx(n :S /2Jn) nes m 
(0, 1), we have 



det(/-Kx( , s /2v^)) < e 



-trK X 



(This follows from the inequality 1 — A < e A when A G [0,1].) The trace 
equals 

>/2y/ri 

K{x, x) dx. 

Now, the usual asymptotics for the Hermite polynomials for "small" x hold 
uniformly for x = 0(n £ ) if e < 1/6 ([1], (19.13.1) and (19.9.5)), so they cer- 
tainly hold for x = 0{^Jn\ogn). Therefore, when s = 0(T,/nlogn), the above 
integral, which equals 

-K[ ) dx, 



2y/Ti V2 v / n'2 v / ra 
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also equals 

1 r s i 
vr Jo 



f s f^ sin2x , \ ^ . 



for some 5. This gives (5.6) and concludes the proof of (5.4). 

For the highest curve, we use the top scaling in trace norm of the m = 1 
kernel, which is just Hn(x,y), as described in Section 5.3. We find, for fixed 
r, that X n (a random variable independent of r) defined by 



V2r(l-r) 2 2 /3nV6 

has the limit distribution F2, the limiting distribution of the largest eigen- 
value in GUE [23]. Since 

we see that as n — > 00, 

(5-7) E (^) = ^^+^6+o(-- 1/6 ) ; 



where 



Here, 



00 

s dF a (a) = -1.771086807411601 . . . 

00 



as computed to high precision by M. Prahofer, so 

cj/~ -0.619623767170.... 

Evaluating (5.7) for n = 5,...,9 gives 2.00981, 2.26104, 2.49069, 2.70345, 
2.90254, respectively. These numbers should be compared with the entries 
in the third column of Table 1. 



Acknowledgments. The authors thank Chris Cosgrove for his invaluable 
assistance regarding the differential equation (1.2) and Neil O'Connell for 
elucidation of the other Bessel process. The first author thanks Anne Boutet 
de Monvel for her invitation to visit, and her kind hospitality at, the Uni- 
versity of Paris 7, where part of this work was done. 



NONINTERSECTING BROWNIAN EXCURSIONS 27 
REFERENCES 

[1] Abramowitz, S. and Stegun, L. A., eds. (1964). Handbook of Mathematical Func- 
tions with Formulas, Graphs, and Mathematical Tables. U.S. Government Print- 
ing Office, Washington, DC. 

[2] Aptekarev, A. L, Bleher, P. M. and Kuijlaars, A. B. J. (2005). Large n limit 
of Gaussian random matrices with external source. II. Comm. Math. Phys. 259 
367-389. MR2172687 

[3] Brezin, E. and Hikami, S. (1997). Spectral form factor in a random matrix theory. 

Phys. Rev. E (3) 55 4067-4083. MR1449379 
[4] COSGROVE, C. M. (2006). Letters to authors, June 9-10. 

[5] Cox, D., Little, J. and O'Shea, D. (1997). Ideals, Varieties, and Algorithms: An 

Introduction to Computational Algebraic Geometry and Commutative Algebra, 

2nd ed. Springer, New York. MR1189133 
[6] Daems, E. and Kuijlaars, A. B. J. (2005). Multiple orthogonal polynomials 

of mixed type and non-intersecting Brownian motions. Preprint. Available at 

arXiv:math.CA/0511470. 
[7] Erdelyi, A., ed. (1953). Higher Transcendental Functions 2. McGraw-Hill, New 

York. MR0058756 

[8] Eynard, B. and Mehta, M. L. (1998). Matrices coupled in a chain. I. Eigenvalue 

correlations. J. Phys. A 31 4449-4456. MR1628667 
[9] Fisher, M. E. (1984). Walks, walls, wetting, and melting. J. Statist. Phys. 34 667- 

729. MR0751710 

[10] Flajolet, P. and Louchard, G. (2001). Analytic variations on the Airy distribu- 
tion. Algorithmica 31 361-377. MR1855255 

[11] Ito, K. and McKean, H. P. (1974). Diffusion Processes and their Sample Paths. 
Springer, Berlin. MR0345224 

[12] Jimbo, M. and MlWA, T. (1981). Monodromy preserving deformations of linear or- 
dinary differential equations with rational coefficients. II. Phys. D 2 407-448. 
MR0625446 

[13] Johansson, K. (2003). Discrete polynuclear growth and determinantal processes. 
Comm. Math. Phys. 242 277-329. MR2018275 

[14] Gillet, K. (2003). Asymptotic behaviour of watermelons. Preprint. Available at 
arXiv:math.PR/0307204. 

[15] Karlin, S. and McGregor, J. (1959). Coincidence probabilities. Pacific J. Math. 
9 1141-1164. MR0114248 

[16] Katori, M., Tanemura, LL, Nagao, T. and Komatsuda, N. (2003). Vicious walks 
with a wall, noncolliding meanders, and chiral and Bogoliubov-de Gennes ran- 
dom matrices. Phys. Rev. E 68 021112. 

[17] Krattenthaler, C, Guttmann, A. J. and Viennot, X. G. (2000). Vicious walk- 
ers, friendly walkers and Young tableaux: II. With a wall. J. Phys. A 33 8835- 
8866. MR1801472 

[18] Louchard, G. (1984). Kac's formula, Levy's local time and Brownian excursion. ,7. 

Appl. Probab. 21 479-499. MR0752014 
[19] Majumdar, S. N. and Comtet, A. (2005). Airy distribution function: From the area 

under a Brownian excursion to the maximal height of fluctuating interfaces. J. 

Statist. Phys. 119 777-826. MR2151223 
[20] Prahofer, M. and Spohn, H. (2002). Scale invariance of the PNG droplet and the 

Airy process. J. Statist. Phys. 108 1071-1106. MR1933446 
[21] Revuz, D. and Yor, M. (1999). Continuous Martingales and Brownian Motion, 3rd 

ed. Springer, Berlin. MR1725357 



28 



C. A. TRACY AND H. WIDOM 



[22] TAKACS, L. (1991). A Bernoulli excursion and its various applications. Adv. in Appl. 

Probab. 23 557-585. MR1 122875 
[23] Tracy, C. A. and Widom, H. (1994). Level-spacing distributions and the Airy 

kernel. Comm. Math. Phys. 159 151-174. MR1257246 
[24] Tracy, C. A. and Widom, H. (1994). Level spacing distributions and the Bessel 

kernel. Comm. Math. Phys. 161 289-309. MR1266485 
[25] Tracy, C. A. and Widom, H. (1994). Fredholm determinants, differential equations 

and matrix models. Comm. Math. Phys. 163 33-72. MR1277933 
[26] Tracy, C. A. and Widom, H. (2004). Differential equations for Dyson processes. 

Comm. Math. Phys. 252 7-41. MR2103903 
[27] Tracy, C. A. and Widom, H. (2006). The Pearcey process. Comm. Math. Phys. 

263 381-400. MR2207649 



Department of Mathematics 
University of California 
Davis, California 95616 
USA 

E-MAIL: tracy@math.ucdavis.edu 



Department of Mathematics 
University of California 
Santa Cruz, California 95064 
USA 

E-MAIL: widom@ucsc.edu 



